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ABSTRACT 

It is well-known that the peak brightness of the Type la supernovae calibrated with Cepheid 
distances can be used to determine the Hubble constant. The Cepheid distances to host 
galaxies of the calibrating supernovae are usually obtained using the period-luminosity (PL) 
relation derived from Large Magellanic Cloud (LMC) Cepheids. However recent empirical 
studies provide evidence that the LMC PL relation is not linear. In this Letter we deter- 
mine the Hubble constant using both the linear and non-linear LMC Cepheid PL relations 
as calibrating relations to four galaxies that hosted Type la supernovae. Our results sug- 
gest that the obtained values of the Hubble constant are similar. However a typical error 
of ~ O.OSmag. has to be added (in quadrature) to the systematic error for the Hubble con- 
stant when the linear LMC PL relation is used, assuming that the LMC PL relation is in- 
deed non-linear. This is important in minimizing the total error of the Hubble constant in 
the era of precision cosmology. The Hubble constants calibrated from the linear and non- 
linear LMC PL relation are = 74.92 ± 2.28(random) ± 5.06(systematic) km s"^ Mpc"^ and 
Hq = 74.37 ± 2.27(random) ± 4.92(systematic) km s~^ Mpc~^, respectively. Hubble constants 
calculated using the Galactic PL relations are also briefly discussed and presented in the last 
section of this Letter. 

Subject headings: distance scale — Cepheids — galaxies: distances and redshifts 



1. Introduction 

Due to their intrinsic brightness at maximum, it is well-known that the Type la supernovae (SNIa) 
can be used to obtain the Hubble constant (i?o)- Furthermore, the peak brightness of the Type la SN is 
regarded as a standard candle after correcting for their light curve shape (for example: Ami5(i?), s-factor, 
MLCS and CMAGIC in Phillips 1993; Perlmutter et al. 1997; Riess, Press & Kirshner 1996; Wang et al. 
2003, respectively). Nevertheless, the peak brightness of the supernovae after the light curve shape (and 
extinction) corrections need to be calibrated with nearby samples before they can be applied to derive the 
Hubble constant. 

A common way to calibrate the peak brightness of SNIa is by using the period-luminosity (PL) relation 
from Cepheid variables in the galaxies that host the supernova (see, e.g., Gibson et al. 2000; Saha et al. 
2001; Riess et al. 2005, and the reference therein). The current most widely used PL relation is derived 
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from Large Magellanic Cloud (LMC) Cepheids. For a long time, the LMC PL relation has been regarded 
as linear in log(P), where P (in days) is the pulsation period of the Cepheids. However, recent empirical 
studies have implied that the LMC PL relation is not linear: the LMC PL relation can be broken into two 
relations, for short (log[P] < 1.0) and long period LMC Cepheids, respectively (Tammann & Reindl 2002; 
Kanbur & Ngcow 2004; Sandagc, Tammann & Rcindl 2004; Ngcow ct al. 2005; Kanbur & Ngcow 2006). 
Various rigorous statistical tests have been performed and the results strongly suggest that this non-linearity 
is real and not due to other factors such as extinction errors or a small number of long period Cepheids^. 
Therefore, it is of great interest to examine how the non-linear LMC PL relation affects the calibration of 
the Hubble constant. 

Currently there arc two studies that deal with this problem: Ngcow & Kanbur (2005, hereafter NK05) 
and Riess et al. (2005, hereafter R05). In the former study, the authors examined the linearity of the LMC 
Wesenheit function, a linear combination of PL and PC relations. The Wesenheit function is frequently 
applied to derive Cepheids distances because it is reddening free (see. e.g. Preedman et al. 2001; Saha et 
al. 2001; Kanbur ct al. 2003; Leonard ct al. 2003). NK05 found that the Wesenheit function for the LMC 
Cepheids is linear because the non-linearity of the LMC PL and PC relations almost cancel out. NK05 
also suggested that the effect of a non-linear LMC PL relation in distance scale applications is minimal (at 
~ 0.03mag. level). However, the authors did not go a step further to compare the Hubble constant calibrated 
from the linear and non-linear LMC PL relation and show that this is indeed the case. In contrast, R05 used 
the long period part of the non-linear LMC PL relation (their "OGLE-l-10" PL relation) to calibrate the peak 
brightness of the SNIa and hence derive the Hubble constant. But their study lacks a detailed comparison 
of the effect of linear vs. non-linear PL relation in distance scale applications. Therefore the main purpose 
of this Letter is to bridge the gap between these two studies by comparing the Hubble constants calibrated 
from the linear and non-linear LMC PL relation. 



2. Data &: Analysis 

Following the prescription given in R05, the Hubble constant from SNIa can be obtained with the 
following equation: 

log Hq = Q.2Ml{tmax) + 5 + aA, (1) 

where Ml(tmax) is the extinction corrected absolute magnitude (in bandpass A) at peak brightness, and 
a\ = log(c2:) — Q.2m\{tmax) is the distance- scale free intercept parameter determined from the "distant" 
supernovae that are located well within the Hubble flow (see, e.g., Jha et al. 1999; Reindl et al. 2005). Here, 
we adopt the same value of av = 0.697 ±0.005 as in R05. This value is determined from 38 SNIa in the Gold 
Sample of Riess et al. (2004). Therefore, once the value of My{tmax) is calibrated with a Cepheid distance, 
the Hubble constant can be obtained in a straight forward manner. 

R05 also listed 4 "ideal" SNIa for the purpose of calibrating the My{tmax)'- they are SN 1994ae (in 
NGC 3370), SN 1998aq (in NGC 3982), SN 1981B (in NGC 4536) and SN 1990N (in NGC 4639). Since 
= I^Cephfij where my{tmax) has been corrected for extinction and light curve shape 



Furthermore, the LMC period-color (PC) relation is also non-linear. 
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of individual supernova, and because we are only interested in the changes of /icep/i,o, this equation can be 
re- written as My{tmax) + I^Cephfi = mvi^max) = a constant (with the constant term uniquely determined 
from the observations of individual supernova). Table 13 of R05 has gives the values of My{tmax) and 
fJ-Cephfi for these four SNIa, hence the change of My{tmax) due to the recalibration of ncephfl is just 
M^{new) = M^{R05) + ncephA^^^) - l^iCephA'^ew). 

To obtain nceph.oinew), we use four sets of LMC PL relations as given in NK05^. Each set of PL 
relations contain both the linear and non-linear (i.e., the long period PL relation) version of the LMC PL 
relation. The Cepheid data for these four galaxies are adopted from the following sources: NGC 3370 from 
R05; NGC 3982 from Steton & Gibson (2001); and NGC 4536 & NGC 4639 from Gibson et al. (2000). As in 
Frcedman et al. (2001), we apply a period cut to Cephcids in NGC 3982, NGC 4536 & NGC 4639 to avoid 
the incompleteness bias at the faint end of the Cepheid PL relation (there is no need to do this for NGC 
3370, see R05). After fitting the PL relations to the data, we obtain the distance modulus via the Wesenheit 
function (see, e.g., the reference in NK05), /xq = ~ 2.45(/Uv — /ij). Metallicity corrections to /^o were done 
in the same manner as in R05 (i.e., using the values in their table 8). Finally, the CTE (charge transfer 
efficiency) correction of -0.07mag. is applied to NGC 4536 & NGC 4639 (Gibson et al. 2000; Freedman 
et al. 2001). No CTE correction is needed for NGC 3370 & NGC 3982. Our results of the ncephfi{new), 
My{new), \ogHo from equation (1) and Hq are summarized in Table 1. 

The random errors that contribute to the Hubble constant include the random error in Cepheid distance 

modulus and the error from supernova light curve fit. Since the random errors in distance moduli, as given 
in Table 1, are almost identical to each other when using either the linear or non- linear PL relation (among 
the four sets of PL relations), we adopt a single value for the error for each galaxy. The random errors from 
supernova light curve fits are given in R05 with values of 0.12mag. for each calibrator. The adopted values of 
the random errors for the four galaxies/calibrators are listed in Table 2. The systematic errors are discussed 
below, and they are summarized in Table 2 as well. These systematic errors are: 

1. Distance to the LMC: the PL relations given in NK05 are based on the LMC distance of 18.50 ± 
O.lOmag., as adopted by Freedman et al. (2001), hence we continue to adopt the conservative value of 
O.lOmag. as the uncertainty in LMC distance. 

2. Linear vs. non-linear PL relation: Since there is growing evidence that the LMC PL relation is non- 
linear, we assume that the non-linear LMC PL relation is the true underlying PL relation. Then there 
is additional ~ O.OSmag. systematic error for the derived distance modulus when using the linear 
version of the LMC PL relation (NK05)^. This error is not applicable to the distance modulus (or the 
Hubble constant) when using the non-linear PL relation. 

3. Hubble flow: R05 determines the error of the Hubble flow from their Gold Sample is about 0.025mag., 
which we adopted this value as well. 

4. HST camera zero-points: The zero-point uncertainties in V- and 7-band for HST cameras (both 
ACS and WFPC2) are O.OSmag., and the total uncertainty in Cepheid distance measurements is 
y(1.45cry)2 + (2.45c7/)2 = 0.086mag. for all four galaxies. 



^Except for KNB05 that we have updated the PL relations in Kanbur & Ngeow (2006), and referred this updated version 
as KN06 in this Letter. 

^This can be seen from, for example, Table 1. For a given set of PL relations, the distance moduli from the linear PL 
relation are systematically closer /further than the distance moduli, by ~ O.Olmag. to ~ O.OSmag., from the non-linear PL 
relation among the four calibrators. 
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Table 1. Results of using the linear vs. non- linear LMC PL relation'' . 



NGC/SN 


Nceph 




Linear LMC PL Relation 

IJ-Ceph,0 M^(tmax) log-ffo 


Ho 


Non-Linear LMC PL Relation 


Ho 
















TR02 PL Relations 








NGC3370/SN1994ac 


64 


32 


.193 ± 0.033 


-19.083 


1 


.880 


75.9 


32.222 ± 0.034 


-19.112 


1.875 


75.0 


NGC3982/SN1998aq 


29 


31 


.634 ±0.064 


-19.164 


1 


.864 


73.1 


31.658 ± 0.064 


-19.188 


1.859 


72.3 


NGC4536/SN1981B 


35 


30, 


.810 ± 0.043 


-19.160 


1 


.865 


73.3 


30.834 ± 0.043 


-19.184 


1.860 


72.4 


NGC4639/SN1990N 


14 


31 


.665 ± 0.077 


-19.025 


1 


.892 


78.0 


31.699 ± 0.076 


-19.059 


1.885 


76.7 














STR04 PL Relations 








NGC3370/SN1994ac 


64 


32 


.171 ± 0.034 


-19.061 


1 


.885 


76.7 


32.149 ± 0.033 


-19.039 


1.889 


77.4 


NGC3982/SN1998aq 


29 


31 


.608 ± 0.064 


-19.138 


1 


.869 


74.0 


31.594 ± 0.063 


-19.124 


1.872 


74.5 


NGC4536/SN1981B 


35 


30, 


.784 ± 0.043 


-19.134 


1 


.870 


74.1 


30.770 ± 0.044 


-19.120 


1.873 


74.6 


NGC4639/SN1990N 


14 


31, 


.648 ± 0.076 


-19.008 


1, 


.895 


78.5 


31.615 ± 0.077 


-18.975 


1.902 


79.8 
















KN04 PL 


Relations 








NGC3370/SN1994ae 


64 


32, 


.189 ±0.033 


-19.079 


1, 


.881 


76.0 


32.180 ± 0.034 


-19.070 


1.883 


76.4 


NGC3982/SN1998aq 


29 


31, 


.627 ± 0.064 


-19.157 


1, 


.866 


73.5 


31.618 ± 0.064 


-19.148 


1.867 


73.6 


NGC4536/SN1981B 


35 


30, 


.803 ± 0.043 


-19.153 


1, 


.866 


73.5 


30.794 ± 0.043 


-19.144 


1.868 


73.8 


NGC4639/SN1990N 


14 


31, 


.663 ± 0.076 


-19.023 


1, 


.892 


78.0 


31.656 ± 0.076 


-19.016 


1.894 


78.3 
















KN06 PL 


Relations 








NGC3370/SN1994ae 


64 


32, 


.195 ± 0.034 


-19.085 


1, 


.880 


75.9 


32.213 ± 0.034 


-19.103 


1.876 


75.2 


NGC3982/SN1998aq 


29 


31, 


.633 ± 0.064 


-19.163 


1, 


.864 


73.1 


31.643 ± 0.065 


-19.173 


1.862 


72.8 


NGC4536/SN1981B 


35 


30, 


.808 ± 0.043 


-19.158 


1, 


.865 


73.3 


30.819 ± 0.042 


-19.169 


1.863 


72.9 


NGC4639/SN1990N 


14 


31, 


.670 ± 0.076 


-19.030 


1, 


.891 


77.8 


31.697 ±0.076 


-19.057 


1.886 


76.9 



^Errors in fio arc random errors only. Unit of Hq is in km s"^ Mpc"\ TR02 = Tammann & Reindl (2002); STR04 = Sandage, 
Tammann & Reindl (2004); KN04 = Kanbur & Ngeow (2004); KN06 = Kanbur & Ngeow (2006, this is tiie updated version of KNB05 in 
NK05). 



Table 2. Error budget for the Hubble constant. 



Source 


NGC3370/SN1994ae 


NGC3982/SN1998aq NGC4536/SN1981B 


NGC4639/SN1990N 






Random Error 






Rl: Cepheid distance 


0.034 


0.065 


0.044 


0.077 


R2: SN light curve fit 


0.12 


0.12 


0.12 


0.12 


Total R 


0.125 


0.136 


0.128 


0.143 






Systematic Error 






SI: LMC distance 


0.10 


0.10 


0.10 


0.10 


S2: Linear vs non-linear PL^ 


0.03 


0.03 


0.03 


0.03 


S3: Metallicity correction 


0.03 


0.03 


0.05 


0.08 


S4: Hubble flow(= 5av) 


0.025 


0.025 


0.025 


0.025 


S5: HST camera zero-point 


0.086 


0.086 


0.086 


0.086 


Total S (linear PL) 


0.141 


0.141 


0.146 


0.159 


Total S (non-linear PL) 


0.138 


0.138 


0.143 


0.156 



^Only applicable when using the linear LMC PL relation. 
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5. Metallicity correction: Uncertainties of the metallicity corrections are adopted from table 8 of R05. 
Note that this uncertainty is considered as a systematic error (Freedman et al. 2001; Leonard et al. 
2003) but not random error. 

The random errors for the four cahbrators are used as the weights when calculating the weighted mean 
of the Hubble constant. This procedure also combines the random errors for individual calibrators into the 
overall random error on the Hubble constant. The overall systematic error in Hubble constant is adopted as 
the straight average of the systematic errors for the four calibrators. The results are shown in Table 3 when 
using the four sets of linear and non-linear PL relation. 

3. Conclusion &; Discussion 

In this Letter we study the effect of linear vs. non-linear LMC PL relations in deriving the Hubble 
constant using SNIa. It can be seen from Table 3 that the Hubble constants obtained from either the linear 
or the non-linear PL relation are consistent with each other (within the total error). The difference in the 
value of the Hubble constant obtained from the two methods is ~ 1.3% or less. Our results arc consistent with 
the finding of NK05 that the non-linear PL relation has an minimal effect, at the ~ 0.03mag. or at ~ 1-2% 
level, on distance scale studies for deriving the distance modulus and/or the Hubble constant. Assume that 
the true underlying LMC PL relation is indeed non-linear, with the latest evidence from Ngeow & Kanbur 
(2006), using the linear PL relation will change the distance modulus and hence the Hubble constant at 
the same ~ 1-2% level as compared to using the non-linear PL relation. Since the four sets of PL relations 
used in Table 3 are not totally independent of each other (they share the same LMC Cepheid data from the 
Optical Gravitational Lcnsing Experiment), the Hubble constants in Table 3 cannot be averaged. Here we 
adopt the Hubble constats from KN06 as our final results. Hence, 

Ho{irom linear PL) = 74.92 ± 2.28 (random) ± 5.06 (systematic) km s"^ Mpc"\ 
ifo(from non - linear PL) = 74.37 ± 2.27 (random) ± 4.92 (systematic) km Mpc~^ 

Note that our Hubble constants are consistent with the Hq Key Project result (Freedman et al. 2001, 72 ± 
8 km Mpc~^) and the recent release of the 3- years WMAP result (Spergel et al. 2006, 73.4t|^ km s'^ Mpc"^). 

It is well-known in the cosmology community that there is a degeneracy between Otot (or 1 — ilk) and 
Hq (or h), as shown in Tegmark et al. (2004). Hence in the era of precision cosmology, it is important to 

Table 3. Hubble constant from various sets of LMC PL relations^ . 



PL sets 


Linear LMC PL 


Non-Linear LMC PL 


% of variation 


TR02 


74.96 ± 2.28 ± 5.07 


74.02 ± 2.26 ± 4.90 


1.3% 


STR04 


75.75 ± 2.31 ± 5.12 


76.47 ± 2.33 ± 5.06 


0.9% 


KN04 


75.13 ±2.29 ±5.08 


75.44 ±2.30 ±4.99 


0.4% 


KN06 


74.92 ± 2.28 ± 5.06 


74.37 ±2.27 ±4.92 


0.7% 



^In unit of km s ^ Mpc ^, first and second errors are the random and 
systematic errors, respectively. 
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minimize errors on the estimation of Hubble constant from observation (to less than few percent level) to 
break this degeneracy. In this Letter we concentrate on studying the contribution to the error only from 
the form of the PL relation (linear vs. non- linear), and find that this could introduce an additional error 
at ~ 1-2% level to the total error. Table 3 suggests that although the Hubble constants and the associated 
random errors arc similar when using the linear vs. non-linear PL relation, the systematic errors arc larger 
when using the linear PL relation if the LMC PL relation is indeed non-linear. Hence it is important to 
eliminate this additional error to improve the measurement of the Hubble constant to within few percent 
level. However there are other systematic errors that contribute to the total error, such as the uncertainty 
of the LMC distance (O.lOmag.) which remains one of the largest systematic error in estimating the Hubble 
constant. Further refinement of all those systematic errors arc clearly desired. 

NK05 has described a way to derive Cepheid distances with non-linear PL relations when the target 
galaxy consist of both short and long period Cepheids. This can be done by using fj, = ^ ^ ^ J2 l^t'^"^* + 

the distance moduli for short and long period Cepheids respectively. 
Most of the HST observed galaxies only contain the long period Cepheids, hence the long period part of the 
non-linear PL relation can be used. In addition to the distance scale studies, the existence of a non-linear 
PL relation in the LMC is very important for stellar pulsation and evolution studies: it is clearly crucial to 
investigate the underlying physics behind non-linear LMC PL relations (see, e.g., Kanbur & Ngeow 2004, 
2006; Ngeow et al. 2005). 

We emphasize that the PL relations used in this Letter are the (linear and non-linear) LMC PL relations 
only. There are recent studies that suggest the Galactic Cepheids follow a different PL relation than then 
LMC Cepheids (Tammann, Reindl & Sandage 2003; Ngeow & Kanbur 2004; Sandage, Tammann & Reindl 
2004), presumably due to the metallicity effects (see, e.g., Gieren et al. 2005, for the opposite point of view). 
Kanbur et al. (2003) compared the Cepheid distances to 25 HST observed galaxies using both the (linear) 
LMC and the Galactic PL relations and found that there is, on average, a negligible (~ O.OOlmag.) difference 
in the distance moduli when appropriate metallicity corrections are applied to the distance moduli from both 
of the LMC and the Galactic PL relations, which is first applied in Kanbur et al. (2003). A similar result is 
also found in a recent paper by Saha et al. (2006) . This may suggest that the use of the Galactic and LMC PL 
relation could have minimal impact on the Hubble constant if the metallicity correction is applied (again, the 
contribution of using the correct vs. incorrect PL relations to the systematic error may be more important in 
reducing the total error on the Hubble constant). To see the effect of using the Galactic PL relation on the 
Hubble constant, we apply the same data and methodology as in Section 2, except that we replace the LMC 
PL relations with the Galactic PL relations. We adopt the recent Galactic PL relations from Ngeow & Kanbur 
(2004) and the updated version in Sandage, Tammann & Reindl (2004). The resulting Cepheid distances 
(with metallicity correction) , the My (new) , log Hq and Hq for the four calibrators are summarized in Table 
4, with the layout similar to Table 1. The error budget for using the Galactic PL relation is very similar to the 
case of using the LMC PL relation (i.e.. Table 2), except that there is no 0.03mag. systematic error from the 
linear vs. non-linear PL relation. Furthermore, we adopt a conservative error of O.lOmag. for the zero-point of 
the Galactic PL relation (Ngeow & Kanbur 2004; Sandage, Tammann & Reindl 2004; Saha et al. 2006), which 
includes various uncertainties from the distance measurements to individual Galactic Cepheids (for example, 
the open cluster fitting method, the infrared surface brightness method and the parallax measurements). 
The resulting Hubble constants are Hq — 70.91 ± 2.16(random) ± 4.69(systematic) km s~^ Mpc""'^ and 
H„ = 69.60 ± 2.12(random) ± 4.60(systcmatic) km s^^ Mpc"^ with Ngeow & Kanbur (2004) Galactic PL 
relations and Sandage, Tammann & Reindl (2004) Galactic PL relations, respectively. The values of the 
Hubble constant from the Galactic PL relation are lower than those obtained from the LMC PL relations, 
however they are still consistent with each others within the la of the total errors. 
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The apparent discrepancy between our results and the results presented in Kanbur et al. (2003) is due to 
the number of galaxies (25 vs. 4) available in both studies. Kanbur et al. (2003) found that if the average of 
log(P) from the Cepheids in a given galaxy is close to ~ 1.4, then the difference in the distance moduli from 
the LMC and the Galactic PL relations will become negligible after the metallicity correction (see equation 
[5] & [6] in Kanbur et al. 2003). This is also true for an ensemble of galaxies. The average log(P) for all 
25 galaxies in Kanbur et al. (2003) is indeed ~ 1.4 (see their table 14), therefore an average of ~ O.OOlmag. 
difference in Cepheid distances is obtained in their paper. For each of the 4 calibrating galaxies in this Letter, 
the average log(P) is greater than 1.46, this makes the Cepheid distances from the Galactic PL relation to 
be further. Hence the resulted Hubble constant obtained with the Galactic PL relations will be lower than 
the Hubble constant obtained with the LMC PL relations. 

The authors would like to thank N. Suntzeff, R. Kennicutt and the anonymous referee for useful sug- 
gestions. SMK acknowledges the support from HST-AR-10673.04-A. 
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